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The main objective of this study nas been to 
determine whether the principles of linear wave super- 
position could be combined with model tests usine wave 
survey techniques to obtain a practical method of 
derivinz “influence diagsrans", which relate the influ- 
ence of small hull form variations on wave resistance. 
Tests were run in the Jebd towine tank usine @ 5.5 ft 
vestroyeretype form with protuberances added at a number 
of locations over a ranzve of speeds. Results showed 
attempts to apply linear superposition for this purpose 
to be larzely unsuccessful. “his result implies serious 
limitations of the linearized theory of wave resistance 
and its application to develop optimum hull forms, 
Yowever, several of tne results point to the conclusio-s 
that, With sone modification, linear Superposition may 
be successfully applied toward the development of 


influence diavrans. 
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INTRODUCTION 

William Froude's investigation of ship resistance of 
nearly one hundred years azo inferred that waves produced 
by a ship moving through the water are a measure of the 
energy imparted to the fluid. It is only in the past ten 
years, however, that this question has received the close 
attention of hydrodynamicists from an experimental point 
of view. “uch of this attention has been focused on the 
theoretical aspects of the wave resistance as derived from 
the wave pattern of a ship, and experimental verification 
by model tests. Some of the most promising experimental 
techniques developed to date have been assessed by Zecers, 
Sharma, and ward oe with very favorable results. There- 
fore, it is the authors’ contention that the state of the 
art has now been advanced to the point where practical 
applications of wave resistance survey techniques may 
be made. 

An important assumption made in linearized wave theory 
is the purported ability to superimpose the individual wave 
patterns of independent hull form variations to obtain a 
compound wave pattern equivalent to the pattern which would 
have been obtained had the hull and form variations been 


towed simultaneously. 





1. Numbers in brackets indicate references at the end of 


this paper. 
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This paper will report the results of an exhaustive 
investigation into this linear superposition theory, and of 
attempts to develop it into a practical means of predicting 
the wave resistance trends of Destroyer-type hulls. ‘To 
this end the following technique is introduced: 

A Destroyer-type model will be tested to determine its 
bare hull wave slope record by means of the longitudinal 
cut, Wave survey method. A protuberance will then be 
attached to the model and the tests replicated. *rom these 
records the wave resistance of the model may be determined. 
Also, assuming linear superposition of the wave records, 
the wave record of the protuberance alone may be determined 
by: 


CS anhiii 1 teed ls Ymoael 
protuberance bare hull 


That is, the wave record of the protuberance will be assumed 
to be unaffected by its location on the model and to give 
the same wave slope record when attached to the model as 
it would if it were alone ina free stream. This hypothesis 
is investigated in depth in a later section of this paper. 
Tne protuberance wave record will be shifted along the 
model length and added back into the bare hull record. The 
resultant wave record obtained is, by the linear superposi- 
tion assumption, equivalent to the wave record that would 
have been obtained had the model been towed with the 
protuberance at that shifted location on the model. An 
influence diagram, as suggested by Horner (2), can be 


developed by considering the protuberance as an increase of 


Ze 





- 


the sectional area at the point of attachment. An 





influence diagram indicates the change in wave resistance 
per unit change of displacement at the different locations 
on the hull surface, and therefore shows at what points 

a given form should be modified in order to reduoe its 
wave resistance. 

Influence diagrams developed completely from experi- 
mental testing require an exorbitant amount of data re- 
cording and analysis. However, if linear superposition 
methods are successful in this paper, diagram construction 
effort will be greatly reduced, and an extremely useful 


tool will be made readily available to the designer. 
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Introduction 

The particular wave measuring technique to be employed 
in this paper will be the longitudinal cut, wave slope 
method (3), which is selected for the following reasons: 
(a) a longitudinal cut through a wave pattern is readily 
obtained by towing a model past a stationary probe, (b) the 
wave slope record terminates much more rapidly than do 
height records, (c) it is a method recommended by fgeers, 
Sharma, and Ward (1,3), and (da) all of the necessary 
equipment was available at Webb Institute of Naval 
Architecture. 


Basic ‘quations and Definitions 


The wave theory which forms the basis for the lateral 
wave Slope method of survey is developed by Sharma (4) via 
enersy flow considerations across a boundary or continuum. 
The final expression relating the wave resistance to its 


respective wave pattern is given by: 


2 Ceres un, Ben 
con 5 J (Sy *Cy ) wile EB) (33 


-@% 
where Sy and. Cy are Fourier sine and cosine transforms 
of the free wave slope profile. That is, 
ae [_ fy SIN (wx) dx (2) 
oo 
yas . ry COS (WX) 3% (3) 
~ OC) 
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wnere fy is the wave slope corresponding to any given 
value of x, u and w are the longitudinal and transverse 
circular wave numbers respectively, and kp is a velocity 
parameter equal to 2 . To fully understand the steps 
required in arriving at a usable solution of equation (1), 


Simple wave theory will be reviewed. 


¥ UJ) 
0 
A x 
= p— 
PuGeee I 


A representative wave of wave length \. (at point A) 
in the wave pattern of a ship advancing with velocity V,, 
is sketched in Figure 1. From the axes! orientations it 
can be seen that © is the angle between® , the resultant 
coordinate of the wave, and the ship velocity. Two 
conditions must be satisfied: namely, that the wave system 
remain stationary with respect to the ship, and that the 
ship be in deep water. Therefore, wave speed is equal to 

VEY iT . With k termed the wave Hiaheezie/ he and 

letting the subscript zero denote conditions at the ship 


centerline, we have: 


Se 





=Vm COS@ JAN ar = /F/R (ny 
= Vom VE! of 21 =“) A, (5) 


and = A im ® COS% - £sec%e (6) 


From Figure 1 it is seen that 


(= XCOSO+YSING 


a 


“Multiplying equation (6) by equation (7), 
kw =f, (x COSe SEC*O 4+Y¥SINO SEC70) 


KW 
RO 


Roh SECO +k SECO FANG 
wWX +Uy (2) 
where Wek, SECO ; Usk, SECO CANO (9) 


squaring and adding equations (9), 


U= =W\/ W (10) 
W =f nf oUF yig +U (11) 


The circular wave numbers u and w have thus been expressed 


il 


in terms of known quantities and each other, 
A solution of equation (1) is outlined by Ward (2) 


as follows. The wave numbers can be non-dimensionalized 
by dividing each by k.. ° ave 
of  W /_k, = SECO 
Ui = U/ fp = SECOLANO 


Equation (1) then ue ue? to 


d ( 
Kw = af Ja | _ a ee (2- 1/,,'2 
(12) 


Since k, ae 2 272 du : 


2. The primes will be omitted throuzshout the remainder of the 


text. 





Sy is determined in a straightforward manner where x 
varies from zero at some point ahead of the wave to a point, 


N, where the F ry 'S die out. 


(3 SIN 
23 La { _ (wh, x) dx ai 


where each Py corresponds to the wave recorded slope tape 
record. This integration must be carried out for each 
value of u. A plot of SC, (1.e., the wave amplitude) 
versus u, W, or@ gives a measure of the generated energy 
of the wave system. Similar wave patterns will display 
Similar wave amplitude spectra. It can be anticipated 

that two wave records possessing seemingly larce discrepancies 
may have very similar wave amplitude spectra, because of 
the smoothing effect of the integration process. A further 
discussion of these spectra will be conducted in a later 
section. 


Equation (12) can be non-dimensionalized to compute 


ae lo Cay PVm* S » Cyrare (peep en (14) 


The basic equation is now in a final form convenient for 


the wave resistance coefficient, Cw, since 


adaptation of a digital computer program solution. The 
programming details and assumptions will be described in 


Appendix B. 


Te 








Characteristics of Model and Frotuberances 

Stevens model number 393 was modified for use in this 
project by removing the bilge keels and adding turbulence 
stimulators. Figure 2 is the model body plan and Table I 
lists the model and ship particulars. The two protuber- 
ances used are shown in “igure 3. Both are bodies of revo- 


lution and, consequently, all sections are circular. 





TABLE I SHIP = MODEL CHARACTERISTICS 
nr 

Length of load waterline 342.5 ft 
Length between perpendiculars 341.0 ft 
Beam 35.9 ft 
Draft i247 ae 
Displacement - tons salt water 2,296.0 
Wetted surface 14,160.0 ft* 
Cy Ore 
Displacement-Leneth ratio 57 
Bean-Draft ratio 2eas 





TABLE I SHIP-MODEL CHARACTERISTICS (con. 





MODEL 
Lambda ratio 62.0 
Length load waterline 5.524 rt 
Length between perpendiculars Fase “Ft 
Displacement-pounds fresh water 21.04 
Wetted surface 3.682 ft’ 
Measured freeboard aft focsle break 2.17 in 
Trin none 
Number of stations 10 


Exploratory Testing 


In order to proceed with the proposed objectives, it 
was necessary to determine approximate towing speeds as well 
as protuberance sizes which would generate significant wave 
Slope difference signals, Fractical limitations confined 
the size of the protuberances between the design water line 
and the keel. Moreover, selection of a suitable hull location 
was necessary for comparison of the difference signals and 
the shifting technique which would be developed. Regarding 
towing speed, it was felt that the selection should depict 
the typical operational speeds of a Destroyer, since design 
techniques (influence diagrams) were of concern. A three 
Wire resistance gage and Brush recording device were used to 
obtain the wave slope records. Equipment details are described 


in Appendix A. 
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Preliminary model runs were made with two hemispherically 
shaped protuberances attached at station 1 of a 10-station 
model. It was assumed that this forward location would give 
the best repreSentation of a bulb generated wave since it had 
the least hull separation of the protuberances. The results 
of these equipment-model familiarization tests tentatively 
verified that the speeds and the selected protuberance size 
would give sufficient difference signals necessary for 
meaningful computer analysis. Consequently, qualitative testinz 


could proceed. 


iesistanoce Tests 

A standard model resistance test was made to obtain 
sufficient data to develop model residuary and total resistance 
curves. The computed values of C, could therefore be compared 
with C, and C, to determine the relative magnitude of the wave 
resistance at various speeds. It could also serve as a check 
to insure that the computed wave resistance constitutes a 
reasonable percentage of residuary teuietands 

From the results of the aforementioned investizations, 
it was decided to tow off-center at Froude numbers of .242, 
.363, and .484 (corresponding to ship speeds of 15, 22.5, 
and 30 knots respectively), with the larger of the two pro- 


tuberances attached successively at stations 1, 3, and 5, and 
#3 The authors know of no previous wave resistance testing 
of a Destroyer-type model. Therefore, the expected rance 


of values is uncertain. 


l2e 


_— 


the smaller attached only at station 3. To provide addi- 
tional data for translation analysis techniques, later 

tests were conducted with the large protuberance at stations 
2and 3. The protuberances would be investizsated at a 
common waterline for comparative purposes, "igure 4 shows 
the model being towed at Froude numbers of 0.242 and 0,424 
with the large protuberance attached at station 2, 

Any errors introduced when comparing the records at 
each configuration, and the sensitivity of the testing 
equipment must be minimized throughout testing. Additionally, 
the influence on wave resistance of the transverse distance 
between the gage and the model must be consistent. A 
detailed error analysis will be presented in Appendix D. 
Suffice it to say, the followins procedures were employed 
to reduce areas of possible error: 

1. A calibration factor relating the recorded valves 
in mm to a slope in radians is essential for each tape 
record. Since the equipment sain settings tend to drift 
during the testing, calibration records were obtained 
prior to and immediately after each testing configuration, 
The gage was physically rotated throush known angles and 
the ensuing deflections were recorded on Brush recorder 
tape. The values were plotted to determine an averace 
calibration factor for each specific configuration. It 


was noted that these values were a definite function of 
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the equipment settings and that they greatly affeot the 
final wave resistance coefficients. 

2. The employment of micro-switches established a 
common wave slope reference point and a rougn verification 
of the model speed. One micro-switch was positioned on 
the tow rail at the point where the bow at design water 
line passed the wave slope probe, the other was positioned 
at a point near the end of the run. The switches produced 
a set of marks on each Brush recorder tape record, and by 
aligning these visually, the records at a given Froude 
number could be quickly compared for consistency. 

3. Bare hull model tests were conducted at the 
middle speed (Fr = .363) with the distance between the 
model and wave probe being varied between 18 and 24 inches, 
in order to determine a range of off-set distance where 
Cy, remained essentially constant and near the maximum. 

The midpoint of this range was utilized throushout all 
towing configurations. Detailed results are shown in the 
data analysis section. 

Three wave slope records were obtained at each Froude 
number for each of seven model-protuberance configurations, 
each speed varying no more than * .01 ft/sec over a 35-foot 
measured run. Figure 5 shows a sample wave record of the 
bare hull model at Froude number .363 and Figure 6 shows 
& sample wave record with the large protuberance located 
at station 3 at Froude number .363. It is worthy of note 


that all testing (total of 63 acceptable reoords) required 















i. ee 


Pee Wave SLOPE RECORD 
moo Nae ULL M00 EL 
Pim” Once 


- = ae Tn oe ff. a , ar 
pis PEED TE 
MV Pi pf] 

















LEUPET PETS TPT 
gee it a “= : 


“wns 







SE PAPE 





































ee ea PEPE PERS 

‘a ee eer 

: nae ae Hert aa DC ae eal 
oeSEaeSe ee aa. 

f ACE EER fet See an rt 
Gaee MAAY ULES SERRE ae 


te. 
when» 
a: j 
=, 
La 


eae EE et SU 
a a | op et Pe Ee ee Sf Ht coal 
a ge EE ee garrett nana 
#E SEE a rol ae 
= Ae 
EE en ba Mic ciao Sitte VEER 
=. = 


aid 

























FIGURE 6 
SAMPLE WAVE SLOPE RECORD 
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only an eleven-hour period. This time included preparing 
the model for each confisuration, reballasting, trimmine, 
and calibration of the wave slope measuring probe, as well 


as conducting the series of off-center runs, 
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AETHOD Of DATA ANALYSIS 

Since there have been no known previous wave resist-~- 
ance tests on Destroyer-type models, a method of reading 
and processing the data had to be developed before any 
detailed analyses, comparisons, or protuberance shifting 
could be initiated. Three representative sets of data were 
used to standardize the data reading and processing pro- 
cedures. Data used were for the bare hull confisuration 
with the larze bulb attached at station 1 at a Froude 
number of 0.242, The following discussion pertains to 
these three sets of data, unless otherwise specified. 

The records obtained from the model tests had to be 
magnified before they could be accurately read, Several 
methods were tried. It was found that the easiest and 
most reliable method was to project the wave slope records 
on the wall next to the computer console and to read the 
data directly onto computer tape. the accuracy of the 
readings thus obtained is estimated to be * 0.2 mm. 

Two alternatives were available for making use of 
the three data tapes obtained for each speed and protuberance 
location. Data tapes could be averased and a wave resist- 
ance coefficient calculated from the resultine data, or 
a resistance coefficient could be obtained for each data 
tape, and the coefficients averavzed for a final wave 


resistance coefficient. Goth methods were applied and 


ia 





the difference between wave resistance coefficients 

obtained was sufficiently small as to be attributable to 

the rounding-off error introduced in averaging the indi- 
vidual data points. However, there was sufficient difference 
between the wave resistance coefficients obtained from the 
individual tapes to preclude use of a single tape as 
representative of the given configuration. 

It was decided to analyze the wave records on the 
basis of an averaged data tape since this method adapts 
itself more easily to shifting the protuberance and requires 
less computer core time. 

Two reading intervals were considered, 2.5 mm and 
5.0 mm. Data tapes for each configuration were read at 
2.5 mm. intervals. In each case the tape length read was 
600 mm. Wave resistance coefficients were then obtained 
by using (a) all data points obtained (240 points, 2.5 mm. 
interval), and (b) every other point (120 points, 5.0 mm. 
interval). Comparison of the wave resistance coefficients 
for a given configuration and speed showed less than 3° 
deviation, as shown in Table II. It was therefore decided 
to use the 5.0 mm, interval in further analyses to simplify 
the tape reading procedures, This will have no effect on 


the results since they are comparative in nature. 
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Jave Sesistance Coefficients at 2. 2 Intervals 





aot. Protub. Frotub., Length Wo. of C, at Cy at & Dart. 


Size Loca. of Data Inter- 2.5 mm. 5.0 mm. from 
Tape vals Inter=- Inter- 5.0 mm. 
nead Orr u vals vals value 
(Sta.) (mm.) (x107) (x10?) 
0.242 Large 1 240 100 O5.317 O23 -~ 1.9 
0.363 Bare - 240 100 0.648 0.666 ~ 2,7 
Hull 
0.434 Bare 
Hull - 240 100 2.359 2,344 + 0.6 


Theory states that the intercration of wave slope data 
with respect to the transverse circular wave number cover 
the range -co to +o , The basic interval between wave 
numbers used was 0.1 (See Appendix C) and resistance coeffi- 
clents were obtained for lengths of 0 through 125 intervals 
of u. The results showed that beyond 100 intervals the increase 
obtained in the resistance coefficient was less than 1%. 

To determine the maximum number of data points re- 
quired to obtain reliable wave resistance coefficients, two 
procedures Were used. First, the data tapes were analyzed 
for lengths varying between 400 and 600 mm. The results 
did not appear to be approaching a limiting value for any 
of the three froude numbers. Therefore, the tank geometry, 
figure 7, was developed, assuming the bow wave angle to be 


equal to the Kelvin angle (194°), and a maximum tape length 
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to avoid reflection was determined. The data tapes were 

then rew-analyzed. Figure 8 is a graphical representation 

of the final analysis. The final data tape reading lengths 
selected were 800 mm. for a Froude number of 0.242, since 

the wave records die out beyond this tape length, and 700 mn. 
and 550 mm. for Froude numbers of 0.363 and 0.434 respectively, 
since reflection occurs beyond these lengths. 

The data obtained for model tests with varying offset 
was analyzed for a data tape length of 600 mm. and a reading 
interval of 5.0 mm. The results, shown as Figure 8, are 
considered to verify that the offset used in the main 
resistance tests (224 inches) was well within the range 
which gives little or no change of wave resistance coeffi- 


cient with change in offset. 
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DISCUSSION OF RESULTS 

The results were compared in a number of different 
ways. These were: comparison of residuary and total re- 
sistance to wave resistance, linear shifts of the difference 
Signals, translation of the influence diagrams, small 
interval shifts of the difference signals, wave spectra 
analysis, comparison of difference signals, and a comparie- 
son of the effects of protuberance size. Each of these 
and representative curves will be discussed in the 


following section. 


Comparison of Residuary and Total Resistance with Wave 


resistance 

In order to obtain a quantitative comparison of the 
wave resistance component with the total and residual 
components, and thus evaluate the predicted savings in 
resistance in later discussions, “igure 9 has been prepared. 
Table III is also presented to give a percentage comparison 
between the various resistances, These ratios are for the 
most part in keeping with the results of wave resistance 
tests on merchant type hull forms and are thus considered 
as an indication of the validity of the tests performed 


for this report. 











rll ll ro 


4) 


Oxe, 


; 
ann 

pn" :) 

“2G 


| + 


a Oy 
x} 
wi 


VY 
ine 


rn 








a CUVEE BASED 
ON ONLY THREE 
DO'INTS 
Eom 0) 50 20 GO a 20 
MODEL SEED (FFS) 
FIG 9 


en Ce VOOR IN 2 Ole ii 
Bee aCe Nd 
26.6 





TABLE IIT 


Resistance Coefficients versus Froude Number for Bare 


Hull iodel 
Ship ¢ G 
i Speed Cy C.. Cy %. We, 
Kts % % 
x10? x102~ = x10? 
0.242 520 0.204 0.92 5.00 2e,2 a 
0.363 22.5 O02 70) (lave yao SoD i2.6 
0.484 30.0 2.445 3.97 7.53 61.6 32.5 
TABLE IV 


Wave Resistance Coefficient x103 versus Conficuration 


Bulb Bulb Froude Numbers 

Location Size 
(Sta) 0.242 0. 363 0.4384 

Bare Hull - 0.204 0.703 2.244 
1 Large 0.329 0.731 1.956 
2 1 OUZ1e 0.741 2ieec 
3 i Ones 0.744 23623 
4 lt 0.261 0.715 2.473 
5 | Om255 0.641 Zea 
3 small 0.220 0.549 29-520 


Construction of Influence Diagrams 


Table IV gives the computed wave resistance coeffi- 
cients versus Froude number for each of the seven configura- 


tions. Disregarding the small protuberance valves at this 


aT 





- 


time, and dividing each large protuberance coefficient by 
its corresponding bare hull coefficient, influence diagrams 
showing the effect of the increased sectional area at 
various points alongs the hull may be constructed for each 
speed. Figure 10 represents the influence diagrams for 
the model, based completely on experimental values. The 
remainder of this section discusses attempts to replicate 
*“{gure 10 utilizins the linear superposition hypothesis. 
If the techniques applied are successful, the amount of 
model testing and wave record analysis required to obtain 
an influence diagram will be greatly reduced, 

In order to apply the linear superposition hypothesis, 
the PROTUS computer program was modified and adapted to 
the NSRDC 7090 computer. The difference signals were 
shifted one data point at a time with respect to the bare 
hull wave record and the wave resistance coefficient was 
computed for each altered configuration. This is equiva- 
lent to a series of discreet shifts of the protuberance 
along tne model hull. figures 11 and 12 are the non- 
dimensionalized influence diagrams derived in this manner 
for the model with large protuberances attached. The 
experimentally-obtained wave resistance coefficients at 
each station are also noted. These curves show that the 
predicted wave resistance coefficients at a shifted position 


are highly dependent upon the initial position of the 
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protuberance and that the predictions do not always agree 
with the actual wave resistance coefficient. 

At this point the physical events which occur when 
the protuberance is moved longitudinally alonz the model 
during a wave resistance test, which are disregarded in 
the idealized linear superposition theory, should be 
discussed. first, as the probuberance is moved aft along 
the model, the separation between the halves is increased 
and the amount of model interference between the two halves 
is increased, Secondly, the direction of the streamline 
flow approaching the protuberance differs at each location, 
effectively changing the angle of inolidence of the flow. 
Finally, the form of the model varies at each location. 
All of these factors influence the wave pattern cenerated 
by the protuberance and could be the causes of the | 
aifferences obtained in the above influence diarrams. ‘7o 
minimize these effects, only the protuberance at station 
1 will be used in later shifting applications, since the 
difference signal generated by tnis configuration most 
closely approximates the protuberance in a free strean. 

Salvesen (5) found that waves tend to die out in 
phase as well as amplitude as they propagate, and conse- 
quently their phase at a distance from the point of senera- 
tion differs from the phase predicted by wave theory. 
Since it 1s expected that this phenomena is occuring in 
the wave resistance tests conducted for this project, a 


translation was applied to the influence diazram in an 
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attempt to obtain better agreement between the actual and 


predicted wave resistance coefficients. 

By inspection of Figures ll and 12 it was decided 
to use a twenty percent aft translation for a Froude number 
of 0.242 and a thirty percent aft translation for a Froude 
number of 0.363. However, for a Froude number of 0.484, 
{t was obvious that translation of the influence diagram 
would be unrealistic, and therefore, no attempt was made 
to apply the technique to this speed. These translations 
yielded wave resistance coefficients at stations 1, 3 and 
5 which are within a few percent of the experimentally 
obtained values. 

To verify this technique, additional model tests 
were conducted with the large protuberance attached at 
stations 2 and 4 At both Froude numbers these tests 
yielded wave resistance coefficients at station 2 which 
were within six percent of the predicted translated 
values, while the wave resistance coefficients obtained 
at station 4 deviated considerably from the predictions. 
The experimentally obtained points have been shown on 
Figure 13. Table V compares the actual and predicted 


resistance coefficients. 
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Comparison of Actual and Fredicted Wave iesistance 


Coefficients 

fy Large Predicted Actual % Deviation 
Protub or ome of Predicted 
at Sta 3 from Actual Ce 

x102 x10 

0,242 2 O72 31 0.215 6.0 

0.242 ky 0.362 0.261 BOs 7 

0.363 2 Onc 0.741 10.0 

0.363 Ly 0,622 Oa a3. 0 


To check the validity of the linear superposition 
technique for short distances, a comparison was made of 
the actual wave resistance coefficients obtained from 
model tests and the wave resistance coefficients predicted 
by shifting the protuberance one station forward or aft, 
with no translation of the influence diagram, This con- 
parison is shown as Table VI. It is evident that a higher 
froude numbers a shift of the difference signal of one 
Station aft predicts the wave resistance with considerable 
accuracy, while shifting the difference signal forward 


rives no reliable results. 
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TABLE VI 


One Station Shift Comparison for Large Frotuberances 
(a) (2) (3) (4) (5) (6) 


Fr Shift To C, at Cy Pred Diff % Diff 
From (B) at (2) (3) (4) (5) (3) 
0.242 1 2) wie .4373 O45 20.6 
i 3 H  w26l .352 -.091 -34.9 
" 3 2 .218 4380 -.162 miAcieD 
i" 5 *» .264 7295 034 13.0 
0.363 1 2) a. 770 -.029 -3.9 
: 3 4 .715 695 -,020 “2.8 
i" 3 2 .7a 885 -.034 -11.3 
" 5 uy ~715 7/70 -.055 -7.7 
0.484 1 2 2.228 2,184 0.044 2.0 
i" 3 & 2.473 2.510 -0.037 =. 5 
" 3 2 2.428 2,710 -0.482 ~-21.6 
i 5 & 2.473 2.136 0.607 eS 


Data Investication 


At eaoh Froude number under investigation, the wave 
record of the protuberance was obtained by the method 
outlined in the introduction. These wave records or 
difference signals are shown as Figures 14, 15, and 16. 


On each figure the assumed starting location, taken as the 
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protuberance maximum cross section, of the difference 
signal 1s indicated by/ . 

If the original linear superposition theory were 
completely valid, the difference signals for each protub- 
erance location at a given Froude number would be identical 
in phase and amplitude, varying only by their starting 
position. This shift in starting position would be directly 
related to the location of the protuberance on the model. 

As can be seen from the difference sicnal plots, the 
correlation between the curves is not exact. However, 
the following points of agreement should be noted: 

a. All of the difference signals are sinusoidal) in 
nature and are quite similar in character to the wave 
records obtained from the model tests. 

b. For a given Froude number, the difference signals 
exhibit fairly sood phase and amplitude relationships. 

ec. The magnitude of the amplitudes of the signals 
decrease with increased Froude number, indicating that the 
effect of a protuberance upon the wave resistance decreases 
with increased speed. 

d. At all Froude numbers the most severe difference 
Signal occurs for the protuberance at station 3. 

From this discussion it is evident that direct super- 
position may be possible but that additional analysis of 
the data is required. 

The amplitude spectra and “ourler Sine and cosine 


transforms were plotted so that the data could be compared 








after it had undergone the first integration. This process 
sives an indication of the effect of the deviations between 
difference signals upon the wave resistance coefficient 

and allows comparison of the data after it has been smoothed. 
Figures 17, 18 and 19 show the curves obtained at a Froude 
number of 0.363 for the configurations with large protuber- 
ances attached at stations 1, 3 and 5. The remainder of 
the spectra are included in Appendix C. 

The agreement among these spectra is remarkably good, 
particularly at the higher Froude numbers. The resultant 
amplitude spectra and Fourier transform spectra show a 
Similar phase relationship at each Froude number as well 
as definite amplitude agreement. 

In addition to the above,plots, the amplitude spectra 
of the difference signal alone were obtained for the large 
protuberance attached at stations 1, 3, and 5 at all three 
frroude numbers, These curves are shown as Figure 20. The 
nature of these spectra indicatesthat the difference signal 
obtained is essentially equivalent to that of the protub- 
erance in a free stream, since in a free stream the protub- 
erance can be represented by a point source, and thus it 
would give a straight line amplitude spectra. 

Although free wave spectra cannot be directly applied 
to shifting methods, their consistent similarity indicates 
that further investigation of the proposed linear super- 


position technique should be pursued. 
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cize Comparison 

The effect of protuberance size upon the wave re- 
sistance coefficient was analyzed by comparing the results 
obtained frommodel tests with the large and small pro- 
tubverances attached at station 3. The methods of analysis 
were essentially the same methods employed in the 
previous sections and are best explained by the followine: 

a. Difference sicnals, shown as Fisures 21 and 22, 
show no direct relationship. 

be. Wave amplitude and “ourier sine and cosine trans- 
form spectra, shown for Sroude number 0.363 as Figures 
23, 24 and 25 and for Froude numbers 0.242 and 0.434 in 
Appendix C, show relatively sood relationships at the 
hicher ‘roude numbers. 

c. Influence diagrams, *igures 26 and 27, are almost 
identical at *r = 0.434, while there is no agreement be- 
tween the diagrams at the lower Froude numbers. 

ad. The ratio of experimental wave resistance coeffi- 


cients at the three speeds are: 


The Cw Cw 
large small 
OR 232 Cee 
0.363 Vs 
0.484 1.069 


As is evident from tne various curves, the agreement 


between the results is a function of the Froude number 


46. 











with the best asreement occuring at the highest speed. 


This is probably because the relative influence of the 


protuberance decreases with increasing speed. 
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Contrary to the results odtained by <2eers, sharma, and 
ward (1,3) rezarding the feasibility of the lateral cut, wave 
Slope method of determing wave resistance, the authors experi- 
enced larze variations in the computed wave resistance coeffi- 
cients of individual data record sets. .jhen principles of 
statistical probability were introduced (Appendix _), the 
average wave resistance coefficient used was found to have 
only a 50% desree of confidence cf beinzs within 67 of the actual 
value. Therefore, all results in this paper must be viewed 
with a certain desree of uncertainty. 

At the present time, influence diasrams can ve ontained 
only by extensive model testinz, Attempts to accurately 
produce an influence diazram by applyins the pinciple of linear 
Superposition of wave records were met with varyinse devtrees of 
success. The following positive trends were observed: (8; ith 
increasing ‘“roude number the relative effect of the protu-er- 
ance upon the wave resistance decreased, yieldins better 
correlation between the experimental and predicted wave resist- 
ance coefficients. (b) The linear superposition technique of 
developing influence diasrams appears to be valid only over 
small intervals. (c) At the lower Proude nunber, translation 
of the influence diagram was required to ovtain any devrres of 
correlation between experimental and predicted wave resistances. 
Mis technique produced relatively sood azreenent within an 


interval of two stations on the model. 
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vnen two protuberances of different size were tested at 
tne same location, no discernable arreement or uniform trends 
were found at the lower Froude numbers. Only at the ni-shest 
“roude number, where the effects of the protuberance are 
diminished, was any consistency observed, 

Based on the foresoing results, the followinet recomnenda-~ 
tions are made for future work in this field: (a) The statis-~ 
tical nature of the data snould be investigated to determine 
the approximate number of data records required to obtain 
sufficient confidence in the final wave resistance values. 

This appears to be especially true for the type of wave slope 
measurement used in the present investization. (b) cince 

there have been no previous Wave resistance tests on lestrover- 
type models, the percentace relationship between the wave and 
residuary resistances is not known, It is possible that the 

low values obtained in the present tests reflect an error in the 
metnod or instrumentation, *urthner equipment checks in botn 
areas should be made. (c) Curve fitting techniques should be 
applied to the experimentally-obtained Wave resistance coeffi- 
cients to improve the quality of the influence diagsrans, (d)} 
Theoretical investigation into nonelinear wave theory should *e 
conducted to discover, if possible, simple adjustments to the 
linear superposition principle so that acceptable azreenent 
between the experimental and predicted wave resistances and 
between results for different-sized protuberances can be obtained, 
“ne similarities between the difference signals and the wave 
spectra noted earlier encourase this effort. Success would 

allow one to construct a valid influence diazsran for wave resist- 
ance which does not seem possible usins the present linearized 


approach. 
1 26. 
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A suecessful wave slove measurin™ probe based on the 
resistance wire principle was developed in 1965 for use at 
yebb in conjunction with the 2ko excitation from a two-channel 

Srush  J-4622-06 anplifier and ‘ark Ii recorder. The 
probe and the calibration scheme is shown in igure A-1, the 
correspendins ceireult diasram is snown in “icure A-2, and 
pnotorranhs of tne equipment in use are snown as “Leure “#3. 
“he seneme is an adaptation of the wave heisnt measuring prote 
previously developed in connection with this effort (5), 
and eaploys three parallel stainless-steel wires of .004 
inch dianete-, spaced at 3/4 inch intervals between care- 
Fully drilied holes in t[lexislass blocks held by a ~-~shaped 
frane. “he frame is mounted on a slider whicn allows lknowr 
variations of vertical heisht of t 3" at approxinately 
=" intervals and this in turn is on a plate which can be 
rotated in either direction to produce a slope variation of 
i .c0 radians in steps of 0.16 for the slope calibration. 
Me heizht variations are used to calivrate wave height 
probes nounted on the same frame or in the present case to 
check the insensitivity of the wave slope probe to chanres 


in wave elevation. 





4, “his appendix nas been extracted fron " xperimental 
.eternination of ‘ave ‘esistance of a whip ‘odel from 


rT 


lateral jave Jlope ‘easurements" by b. |. hard, cep. 1967 
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As shown in Figsure Ae-2, the measuring circuit consists 
of a Nheatstone bridge with two of the resistances Ry and 
Qn replaced by the equivalent path resistance of the water 
between each set of wires AB and AC. This scheme is known 
in the case of the height probe to give a good linear response 


to changes in wave height over a range of several inches 


when the ratio of resistance: 


is properly chosen; in this case the value of K is approxi- 
mately 20 based on the equivalent D. C. resistance of the 
water elements. In the oase of the slope probe the elements 
are used differentially, as shown. This circuit is more 
symmetric and actually easier to balance than the height- 
probe circuit where an electrical resistor must be used to 
attempt to reproduce the water effect. A combined resistance 
and capacitance balancing circuit is provided within the 
bridge circuit to assist that available in the amplifier. 

As mentioned previously, good linear response has been 
established for this probe in operation over the range of 
* 0.20 radians ample for the usual model wave system and 
relatively good insensitivity with changes in wave height 
over the range of e 14 inches required for the present tests. 
Care is taken to keep the wires clean and at constant tension. 
No difficulty with this arrangement has been noted due to 
variations in surface film, nor does the exact spacing of 


the wires seem to be a critical factor as has been noted by 
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others. This might be a fortuitous result of the relatively 

high exciting frequency of 2KC provided by the Brush amplifier. 
The foregoing statements refer to the characteristics 

noted under static conditions and do not include any effects 

due to dynamic action, such as meniscus csocillation at the 

frequency of the waves being measured, Although Pearlman 

ran a series of such tests on a similar device and found 

no important effects at the frequencies of importance, the 

present device should be subjected to a similar series of 


tests for the sake of completeness. 
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A basic program computerizing Equation (14) is civen 





by Ward (3). An extensively modified version of it was 
adapted for this paper, although the fundamental computa- 
tions have remained intact. In addition to computing a 
wave resistance coefficient, the program (named "PRiOTUB") 
developed for this project can subtract a bare hull record, 
point by point, from a with-protuberance configuration 
record, and compute the difference signal wave resistance 
coefficient. Additionally, it can shift this difference 
record forward or aft along the model, multiply each slope 
value by any desired constant, add the modified sicnal back 
to the bare hull record, and determine the wave resistance 
coefficient of the altered data. 

FROTUB, as listed in Figure Bel, is written in GE 
time-sharing fortran. The input data consists of NZETA> 
wave Slope values which are stored and referenced in the 
SFILE statement. Each data point is divided by the calibra- 
tion factor for that configuration (CAL1 or CAL2) to obtain 
the ZETAY(N) set in radians, which are the proper units 
for processing. 


Equation (14) may be re-written in series form as: 


JU 


gel 2 2 
ceeds By etre aubtledy co 





5. Capitalized words correspond to names as programmed in 


PROTUB. 











The incremental value of \U is given in ief (3) as 


Ax = T/A, bese (B-2) 


where Derr is a hypothetical effective tank width. If 


it is assumed that K barr = 1017, thenA aw =—— = 0.10 c 
in 
N2ETA 
Sy = SY = ZETAY(N) SIN(WA,X) AX 
(3-3) 


Ged 
t 


cy = CY = ¥” ZETAYIN) COS (WR, x) AX 
N>/ 


6, This simplifying assumption may be justified by consider- 
ation of the tank seometry. The actual effective tank 


width considering the model offset equals 11.35 ft. 


we 0.12 
eff : SSS Ve 
4,40 J /? tr 40.5 


Therefore berr is generally equal to or greater than 
the actual value, giving a fu about equal to or smaller 
than required from Equation (8-2). This provides a 


finer integration than is required by Equation (3-2). 
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X = DELX = Viilmodel speed in ft/sec) “ PAP ZR 





(Interval between readings in mm) / VP (paper 
Speed in mm/sec) 
Therefore DELX equals the distance between data 
points in feet 
mG No DEEX, N= 12,395.50 ene A 
wis a funotion of J and is computed from =qua-~ 
tion (11) for each increment of. 
Letting CAYO = oon — factoring DELX out of the expression, 
Dy = z F lalicilimiin (W* CAYO*¥ N* DELX) 
Cy en ZETAY(N) RCOS(WH CAYO#NY DEK) ) 
Pquations (B-4) solve the Fourtler transforms. How- 
ever, the trigonometric identities, 
sin(X +A X) = sinX cosAX + cos sindx 
cos(X #4 X) = cosx cosAhX = sinx sindx 
are utilized in PROTUS to reduce computer core time. 
Letting Kgberr = CAYOB = 1077 


S = WETSR = nodel surface-ft* 


= 2(DELX)® sy? + ov? 
WETSR*CAYOB wr *(2w* 1) (355) 


where the entire éMiresmton ie doubled because only 4 the 
wave pattern is recorded by the probe. Note that the final 
resistance coefficient is the product of a constant and the 
wave amplitude spectra, squared, times a weishtins function, 
1 » which decreases quite rapidly with increase in 


wa? (2w2-2) u or W, 
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APPENDIX C; WAVE SPECTRA PLOTS 

The remaining amplitude and Fourler sine and cosine 
Spectra, not included itn the body of the report, are shown 
in this section. The curves were plotted by the NSRDC 7090 
computer, Due to the limitations of the curve plot sub- 
routines used, the following notation was adopted for identi- 
fication of the curves: 

First letter denotes: 

A Bare hull configuration 


large protuberance at station 


B 
C Large protuberance at station 
D 


large protuberance station 


= Small protuberance station 


First number denotes: 
1 Froude number equals 0.242 
2 Froude number equals 0.363 


3. Froude number equals 0.484 


AVE denotes that the curves are plotted from data averased 


over three model tests. 
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APPENDIX D: ERROR ANALYSIS 





Data Inconsistency Record 
Although the three data tapes retained for each configu- 


ration and speed had an average speed within + 0.01 ft/sec, 
the wave slope records were not identical, and significant 
phase shifts developed toward the middle of most of the 
records. To investigate this source of error four sets 


of data were artibrarily selected. These were: 


set Speed Conficuration 
1 242 Large protuberance at station 1 
2 969 Bare Hull 
3 363 Large Protuberance at station 2 
by 48h Bare Hull 


When individual tapes of each set were analyzed, the 
extreme Wave resistance ooefficient had an average devia- 
tion of almost 10% from the mean. The deviation amoung the 
four sets of data so analyzed are shown in Table DeI. 

TABL® Del DATA DEVIATION 


Cy Of Individual Data Tapes 


Set Tape l Tape 2 Tape 3 Cy Wax Dev from om 

1 23129 3545 .31™ - 3290 yo 

2 »6755 -7780 6016 - 7033 10.7 

3 8404 ~7394 .6790 ©7529 11.6 

bh 2,348 2.382 2.067 2,264 8.8 
Average 9.7 


LOL. 








There is strong reason to believe that the wave records 
are statistical in nature, and that many more data tapes 
are required to attain a suitable degree of confidence in 
the average value, Applying the student 'T' distribution 
for small samples, if Cue and a4 are the mean and variance of 
a sample of size N from an assumed normal distribution, 
nl4,07), where 4{ and 0 * are unknown, then the percent 
deviation of Cay from the mean with a given confidence interval 
may be readily determined. 

Each of the four sets of data individually investigated 
and compared in Table De-I experienced considerable variation. 
Therefore, probable deviations of Cw from the theoretical 
mean were determined for each set at four confidence levels 
and averaged. These steps are tabulated in Figure Del, 
which shows the percent deviation of the accepted Cw from 
the mean versus its respective degree of confidence. For 
comparative purposes, the same values of s and ~ are used 
to calculate the confidence levels were ten samples to be 
taken. It is quite apparent that, from a statistical view- 
point, three runs per set do not produce sufficiently reliable 


values of wave resistance coefficients. 


Calibration Factor 


The Brush recorder and amplifier exhibited some fluctua- 


tion and nonelinearity during testing, A Sample calibration 


factor plot for the bare hull configuration is shown in 


LOZ 










FIGURE D-l CONFIDENCE LEVELS FOR WAVE RESISTANCE COEFFICIENTS 






S/ {3 % DEVIATION FROM MEAN AT % CONFIDENCE 
90 95 


23290 .0143 3255 8.2 ae 18.7 















0/033 .0337 3.92 9.0 14.5 20.7 
27529 .0471 5-20 11.8 18.3 27.0 
2.244 .0830 3202 7.0 10.8 15.9 
9.0 14.2 20.5 
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FPigure De2. Depending on the manner of securing the probe 
after each rotation when calibrating, an error of = } mn. 
at .2 radians was observed. This gives an expected calibra- 
tion error of 7 1.25 mm/rad. 
As the probe is raised or lowered with no rotation, 
a perceptible scale deflection was detected. However, the 
resulting effeot on the wave record tends to be self-con- 
pensating and this factor is mentioned for information only. 
Therefore, since the average calibration factor is 
about 95 mm/rad and is squared in the calculations, per-~ 


centage error 1s approximately (t 1.5/95)°= 23%. 


Wave Record Reading Error 


Because the human element is involved in reading the 
data tapes, a reading error of measurable magnitude is likely 
to be incurred. To investigate this problem, each of the 
tapes in the fourth set of data mentioned above was read 
individually by each of the. authors and separate values of 
Cy computed. Table D-II summarizes the results. From this, 

a gross reading error of . 1 3/4% is anticipated. 


Table DeII Reading Error 


Protuberance located at station 2, Froude number = 0,363 


Reader lst Tape 2nd Tape 3rd Tape Averace 
FIM ~ 3337 = (ann 6732 -7e09 
CAB . 8333 7648 6712 «7500 
AJJ 3492 7317 6926 7448 


Average for all readings = .7409 


Maximum deviation from average = oe = 1,79% 
o 7409 ‘ | 


10h. 





FIGURE D-2 BARE HULL CALIBRATION FACTOR 
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verall Sensitivity of Method 


If a sufficiently larze sample of wave records were 


taken, the percent deviation of C,, at 90% confidence could 


readily be reduced to a few percent. Maximum combined 
calibration and reading errors constitute about 2: + 1 3/4 
or 4% error. Therefore, the longitudinal cut method, when 
conscientiously followed, with an adequate number of records, 
should provide wave resistance values accurate to within 

6% or less. However, when only three samples are used, a 


6% accuracy is available with only a 50% confidence level. 
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